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Abstract 

A consistent theory of the ground state energy and its sphtting due to the pro- 
cess of tunnehng for the Lipkin model is presented. For the functional integral 
in terms of the spin coherent states for the partition function of the model we 
accurately calculate the trivial and the instanton saddle point contributions. 
We show that such calculation has to be perfomed very accurately taking 
into account the discrete nature of the functional integral. Such accurate 
consideration leads to finite corrections to a naive continous consideration. 
We present comparison with numerical calculation of the ground state energy 
and the tunneling splitting and with the results obtained by the quasiclassical 
method and get excellent agreement. 
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I. INTRODUCTION 



Lipkin, Meshkov and Glick |l|] proposed in 1965 an exactly solvable two-level many- 
fermion model which has been used to test various kinds of many-body theories. The model 
Hamiltonian reads 

H = OSo- ^fiSl + Si) (1) 

where / is the coupling constant, uj is the energy difference between the two levels, 

n 

So = ^ (c+mC+m ~ C_^C_m), (2) 
^ m=l 

1 ^ 

^ m=l 

and f2 is the degree of degeneracy of each level. The single particle states are labeled by the 
quantum numbers ±m. The operators 5*0, S± satisfy the commutation relations of the sm(2) 
algebra. 

The transformation associated with the unitary operator e"^"^" leaves the Hamiltonian 
invariant. If x = f^/^ < the mean field groundstate of the model is classified according 
to the trivial representation of the symmetry group. The symmetry is spontaneously broken 
if X > 1. 

A splitting of the levels is then observed, which is analogous to the occurrence of rota- 
tional bands in deformed nuclei. The model is also of great interest in condensed matter 
physics, since it is related to the anysotropic Heisenberg model. In the strong coupling limit, 
X ^ 1, the level splitting vanishes for odd f2. This behaviour is easily understood in the 
framework of Kramers theorem and is related to the well known phenomenon of tunnelling 
suppression for odd spin anysotropic Heisenberg ferromagnets. 

In this paper we formulate a general physical picture of the calculation of the ground 
state energy and its splitting by the instanton method. Such method was applied to the 
anisotropic Heisenberg model for description of the tunneling of the magnetic moment of 
small magnetic particles with large spin 1^-^,^]. Although the physical picture of the spin 
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tunneling was formulated in these papers the quantitative results of the papers are 
not completely correct. 

We show that special care is required when computing the instanton contributions in 
order to take also correctly into account the small amplitude quantum fluctuations. On 
more technical language the functional determinants have to be calculated very accurately 
taking properly into account the operator ordering generating the functional integral for 
the partition sum. We show that it is essential to keep in mind that the fuctional inte- 
gral is generated by the spin coherent state method. Taking accurately into account the 
discrete nature in time of the functional leads to essential corrections to the ground state 
energy and to its tunneling splitting. After taking into account all these contribution our 
result completely coinsides with the one of obtained by the semiclassical method and the 
contradiction existing in literature is resolved. 

The structure of the paper is as foUws. In Sec. |T| we present the functional integral 
for the Lipkin model. In Sec. |T| we review the the simple results for the tunneling in 
the framework of the continuous approximation for the functional integral. In Sec. |V| we 
calculated accurately the contribution of the trivial saddle point to the partition function 
and determine the ground state energy. In Sec. |V| we calculeted the determinant and 
the measure of integration over time in the framework of the canonical fromulation of the 



functional integral. In Sec. VI we calculate accurately the functional determinant for the 



instanton saddle point taking into account the corrections due to discretization. In Sec. |VI1 
we present comparison between the analytical theory and the numerical calculations for the 
Lipkin model and find an excellent agreement when additional contributions discovered in 
this paper are taken into account. 

II. FUNCTIONAL INTEGRAL FOR THE LIPKIN MODEL 



It is convenient to rewrite the Hamiltonian of the model in the rotated reference frame 
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H = H2 + Hi, Hi = CoSy, (3) 
H, = f{Sl - Si) = f{^-Sl - \s{s + 1) - \{Sl + Si)). 

Here we assume that the magnitude of spin s ^ 1. The first term in the Hamiltonian 
(§) describes interaction between nucleons, the second one represents an energy sphtting 
and can be regarded as an external magnetic field. The choice of the coordinate axis is 
determined by the convenience condition for further calculations. 

The Lipkin model can be reduced to the anisotropic Heisenberg model 0. For that we 
add to the Hamiltonian (^) the term f {S^. + Sy + Si — s{s + 1)) = and the Hamiltonian 
(I) takes the form 

H = -fs{s + 1) + f{2Sl + Si) + uSy. (4) 

In the language of the Heisenberg model the magnetic field here is in the direction of the 
second axis y. The parameter A of the paper is 1/2 for our case. We concentrate our 
attention on the accurate treatment of the ground state energy of the model and its tunneling 
splitting when the second term in @ is absent or u; = 0. 

In this special case the Lipkin model is invariant under the time reflection. In the case of 
half-integer spin s this symmetry leads to the two-fold Kramer's degeneracy of the ground 
state PJ^. In the case of integer spin s instead of degeneracy we have the splitting of the 
ground state. This splitting has to be small because for large s the difference between integer 
and half-integer spin has to be small. In fact this splitting is exponentially small. 

To solve this problem we have used the instanton method ^ applied to the functional 
integral for the partition function of our spin system in terms of spin coherent states [Q. 
The coherent states are given by the following formula 

1^) = (1 + |z|2)-^exp(z5+)|0), 5_|0) = 0, (5) 
^,|0) = -s|0), 

where 2; is a complex number. They possess many remarkable properties and with the help 
of them the partition function can be represented in the form f^: 
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Z ^Tr[ex.p(-TH)], 

r ■■■ r Y[ ^ 

i-OO J-OO „_n nil 



(6) 



exp{A{z)). 



'-oo„=o vr(l + |z„|2)2 

Here the interval of the imaginary time T is spht into parts: T = iVA, and in every 
section an integration over z!^ — ReZn and z'^ — Imzn is performed. The action of the system 
A{z) has the form 



N-l 



A{z)^J2[2sHl + z:^,Zn)- 



(7) 



n=0 



2sln{l + \zn\')-AH{z:^„Zn)], 



where the variables Zn satisfy the periodic boundary conditions: zn = Zq or in the continuum 
hmit z{T) — z{Q), and the Hamiltonian 



{Zn+l 


H 


Zn) 


{Zn+1 


z 


n) 



The matrix elements of all essential operators have the form 

(^l^ol^) = - z*z'){l + z*z')Qsi 

{z\S_\z') = 2sz'{l + z*z')Qs, 

{z\S+\z') = 2sz*{l + z*z')Qs, 

{z\Sl\z') = 2s{2s-l)z*'^Qs, 

{z\Sl\z') = 2s{2s-l)z'^Qsi 

{z\Sl\z') = + z*z'f - 2s{2s - l)z*z')Q, 

{z\z') = {l + z*z'fQs, 

= (1 + z*zT~V{'^ + + \Ay- 



(8) 



(9) 



III. NAIVE CLASSICAL PICTURE FOR THE GROUND STATE ENERGY AND 

ITS INSTANTON SPLITTING 



5 



A. Classical action and its minima 



In many cases (but not in all !) one can take the continuum limit when Zn+i is close to 
Zn- In this continuum limit it is convenient to introduce the canonically conjugated variables 
(fi and p = cos 6: 

S — p < ip < 27T, 

z = pexp{i(p), p = — — , . (10) 

-s<p<s 

The action of the system A{z) in the continuum limit is 

Aiip,p) = r[i{p-s)ip-H{p,ip)]dt, (11) 
Jo 

H{p,ip) = ((?/2)(p2(i + cosV) - s^cosV), (12) 

where H{p, (p) is the Hamiltonian of the problem, g = f{2S — 1)1 S, (7 > is the coupling 
constant. This Hamiltonian is unusual because the mass depends on the coordinate y?. The 
second essential peculiarity is the presence in the action of the term is(p which separates 
integer and half-integer spins . In spite of that the ground state energy and the splitting 
can be found. 

One can easily check that the Hamiltonian has two minima at the points ^9 = 0, vr and 
p = 0. These minima are deep for s ^ 1 and the Hamiltonian in the neighborhood of the 
minima has the simple form of a harmonic oscillator 

H{p, ip) = Emin + gp^ + 5'sW2, Emin = (13) 

and the vibration frequency is = \/2gs. One can assume that the energy of the ground 
state 

Eq naive = Emin + 'hujjl. (14) 

This expression for Eq can be obtained if we calculate the contribution of the trivial saddle 
piont to the partition function by the method of 0. 
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B. The instanton contribution to the spHtting of the ground state 



Besides the trivial saddle point contribution to the partition function discussed above 
(0) there exists a set of nontrivial saddle point contributions. These saddle points, in their 
main features, can be discussed in terms of the continuum action (|rT]). The discussion of 
the Gaussian fluctuations around the saddle points has to be done more accurately but the 
result can be expressed in terms of the continuum action only. The instanton saddle points 
for the action ( ]TT| ) are realized at the imaginary momentum pt = isut- The Hamiltonian 
equations of motion which follow from ([Til) are: 

= 5fs(l + cos^ (15) 
ii = gs sin ip cos Lp{\ + u^). 

We are interested in the time periodic solution of Eqs. (p!5|): (y9(0) = v^(T), j9(0) = p{T). 
These periodic solutions can be constructed from the elementary solutions which are named 
instantons and anti-instantons. It is obvious from Fig. 1 that there are two types of 
instantons (A and B) and two types of anti-instantons {A and B) for our problem. 
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FIGURES 
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FIG. 1. Instanton trajectories in the phase cos{ip) — u space. The curve A corresponds to the 
^-instanton and A antiinstanton, the curve B corresponds to the B-instanton and B antiinstanton. 
cos{(p) is along the horizontal axis, the momentum u is along the vertical axis. 

Each of them is labeled by the time at which cos(v9(to)) = and M(to) = =tl- 

The A-instanton starts at the time t — oo from the point of phase space {(p = — tt, u = 
0) at the time t = to arrives at the point {if = 0,u = —1) and at the time t —>■ oo arrives at 
the point {ip = 7t,u = 0). 

The B-instanton starts at the time t ^ —oo from the point of phase space {(f = — vr, u = 
0) at the time t = to arrives at the point {(f = 0,u = 1) and at the time t — > oo arrives at 
the point ((y? = tt, m = 0). 

The A-anti-instanton starts at the time t — > — oo from the point of phase space {ip = 
7i,u = 0) at the time t = to arrives at the point {ip = 0,u = 1) and at the time t ^ oo 
arrives at the point {p = — vr,^ = 0). 

The B-anti-instanton starts at the time t —oo from the point of phase space {p = 
7r,M = 0) at the time t = to arrives at the point {p = 0,u = 1) and at the time t — oo 
arrives at the point {p = —tt, u = 0). 

The analytic form of these solutions can be easily constructed if we take into account 
the energy conservation during the process of motion H{u, p) = Emm- In this case we have 
the following algebraic connection between p and p: 
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u 



1 — COS^ !f 
1 + COS^ ip 



(16) 



Combining formula 
instantons: 



) and (|T^) we get the equation of motion for instantons and anti- 



= ±5'Sa/ 1 — cos^ v?, 



(17) 



where the sign ± corresponds to the two types of instantons. This equation of motion can 
be easily solved and we have 



COS(y9j 



sinh(ro) 



1 



Ut = T- 



/2 + sinh2(ro)' ' ' cosh(ro) 

where tq = uj{t — to). 

The instanton action Aq counting from the energy minimum is 



(18) 



rT 

^0 = QS^ / (1 COS^ ip)dt 

Jo 



(19) 



^1 — COS^ ip) 



dip = 2sln(l + V2). 



^ (1 + cos^ ip) 

The calculation of the instantons and anti-instantons contribution to the partition func- 
tion is based on the idea of an approximate saddle point. The general trajectory for an 
approximate saddle point is shown in Fig. 2. 
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FIG. 2. The scheme of the consecutive fast instanton and antiinstanton transitions at times 
tl, t2, ... ts- The time t is along the horizontal axis, cos((/9) is along the vertical axis. 
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This trajectory starts from the point (f = n, u = and finishes at the same one. It 
represents a sequence of fast instanton and anti-instanton transitions (A and B types). The 
centers of theses instantons and anti-instantons (points where = ±7r/2, u = ±1) are 
situated at the arbitrary (but ordered points) t/. Notice, that for the exact solution we have 
a periodical lattice of instantons and anti-instantons. The contribution of an approximate 
saddle point to the partition function can be presented in the form (0): 

^a,r = L L n ^(^'+1 - tl) ■ (20) 

1=1 



JO 



l[[Ke-^''-''^^dti][det'{Lsad) 



-1/2 



1=1 



The integral over the positions of the instantons ti in Eq. ( pOD represents a sum of contri- 
butions to the partition function over the manifold with almost the same action Aq. The 
factor if" gives the measure of the integration over degenerate minima 0. Its explicit 
form can be derived if we connect the time ti with the amplitude of the zero mode of the 
quadratic form Lins describing of the fluctuations around of the instanton solution. The 
factor exp(— y4o — isipi) gives the exponent contribution to the partition function. The quan- 
tity Sifi is the Berry phase which follows from Eq. ([l^) for the continuum action. This 
phase SLfi is equal to svr for A-instanton and 5-anti-instanton and — svr for S-instanton 
and y4-anti-instanton. 

The last factor in Eq. (pUD, det'{Lsad), represents the determinant of the quadratic form 
describing fluctuations around the approximate saddle point solution excluding eigenvalues 
corresponding to zero mode fluctuations. We understand the determinant as a product of 
the eigenvalues of the corresponding quadratic form. Our aim is to describe the calculation 
of this determinant in more detail. 

(1). One can represent (\.ei'{Lsad) product 

det'(L,,,) = det(Lo) ■ X, X = det'(L,,,)/ det(Lo), (21) 
where det(I/o) is the determinant of the trivial saddle point which we calculate in Sec. |V| 



(^). It is important to notice that the magnitudes of determinants (let'{Lsad) and det(I/o) 
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depend on the method of regularization of the functional integral for the partition function 
or (what is just the same) from the large eigenvalues of the quadratic forms Lgad and Lq. 
But it is obvious (see Fig. 3) that for Tuj/n ^ 1 the large eigenvalues of the quadratic forms 
Lsad and Lq are the same because during almost all time of motion the variables (p{t) and 
u{t) are close to the trivial saddle points where 99 = 0, tt and m = 0. 




FIG. 3. The effective potential energy for the "particle" in tlie field of many instanton tran- 
sitions. The time t is along the horizontal axis the effective potential energy is along the vertical 



axis. 



This means that the quadratic forms Lgad and Lq differ from each other only during 
the time n/uj <^ T. This means that the quantity X can be presented as the ratio of 
determinants in the continuum limit 



X — det' [Lsad cont)/ det(Lo cont)- 



(22) 



(2). For simplification of the problem it is convenient to integrate the action over the 
variable p. For that we spread the limits of integration to ±00 in the expression ( |Tl|) for the 
partition function. Such spreading of the limits of integration can be justified in the case 



5fAs^ > 1, or s>l/Acj, Acj < 1. 



(23) 



Because Au < 0.1 this leads to conclusion that such integration is justified for s ^ 10. 
We want to stress that this condition includes the time quantization step A. Therefore the 
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expression for the partition function takes the form of a functional integral over the variable 
93, and in the leading approximation over 1/s we have 



2ir /.27r ^ 



e 



Ab+A 



d(pt 



n u . : . (24) 



t=Q\j2T{:^g{l + co??ift) 

Ab = -is [ ipt dt, 
Jo 

(1/2)^ (^ip^/g{l + cos^ (ft) - gs^ cos^ (ft^dt, 



A^p 



where Ab is the Berry phase, A^p is the action of the field (p. We want to remind that the 
measure of integration of the functional integral over configuration space depends on the 
time step A. 

It is natural to introduce a new variable V'f in the functional integral (p4D through rela- 
tions: 

di,t = {g{l + cos\t))-^'''d^t, (25) 

/ {g{l + cos''y))-^'^dy. 
Jo 

The explicit form of ipt as function of (pt is an elliptic function but this explicit form is 
unnecessary for further calculations. In terms of ipt the action A^p takes the canonical form 

A^ = -(1/2) ^[^2 _ ^^g2 (26) 
Jo 

where (pt in this formula is a function of ipf One can check that this action has the instantons 
minima for cosip^ according to Eq. (piSf ) and the instanton action Aq according to (|19|). 



Decomposing the action (|26|) around the multi-instanton minima and using (|25|) we get the 
part of the action quadratic in fluctuations in the form 

Al = -(1/2) J [6iJt Lad cent SiJt] dt, (27) 
Lsad cont = ~^ ^ (V't)' 



n^?) = Y(2cosV? - 1)(1 + cosVt")- 
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The function V"{!f^) plays the role of the potential energy for the problem of fluctuations 
around the instantons (see Fig. 3). For the trivial saddle point cos^ ipt = 1, V"{ip^) = uP' and 
the operator Lq cont has the simple form: 

-f'O cont = + ^ ■ (2^) 

The typical function V"{ip^) is shown in Fig. 3. It is equal to cj^ in almost all region 
< t < T, excluding narrow regions of order l/u around the instanton transitions at points 
t = ti where it has wells of depth —uj'^/2. 

(3). At this point we can calculate the ratio of determinants X (|2^) . This calculation 
is based on the small parameter n/uoT ^ 1 and we can prove that the ratio X is equal 
to the product of the ratio of determinants for each region containing an instanton or an 
anti-instanton: 

n 

X = WXi, X; = det'(L;, cont)/ det (Li, cont), (29) 

1=1 

where the operators L;, ins cont and L;, o cont are defined in the region ((t;+i + ti)/2, (t; + 
ti_i)/2) around the position /-instanton ti on the class of functions which satisfy the zero 
boundary conditions at the points ((tj+i + ti)/2 and (ti + t;_i)/2). This statement can 
be justified by the following arguments. Every determinant det{d) is the product of the 
eigenvalues of the operator d on some class of functions. In our case we can consider the 
class of functions with zero boundary conditions at points ((t^+i + ti)/2 and (ti + ti^i)/2). 
In this case all eigenvalues will be nondegenerated. For example, for the periodic boundary 
conditions they will be twice degenerated what is inconvenient. All eigenfunctions of the 
operator d in the region [0, T] can be divided into two classes: those localized close to the 
instanton positions ti with the "bounded" eigenvalues forming almost degenerated groups 
of n eigenvalues and the propagating eigenfunctions with "continuum" eigenvalues almost 
the same as in the case of the trivial saddle point. It is obvious that the contribution 
of the "bounded" eigenvalues to the total determinant can be presented as the product 
of n contributions. The "continuum" eigenvalues can be calculated in the quasiclassical 
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approximation and can also be split into n parts. This argument proves the validity of Eq. 

(4). Because our problem is reduced to the well-known one for the usual instantons 
we can use results of Coleman [0. The quantities Xi of Eq. ( pQ]) can be calculated if we 
construct the normalized eigenfunction x(t) = Nip(t), with = l/^/Ao, of the operator 
Lsad cont (13) with zcro eigenvalue 0. The function x{t) has a form 



X{t) = -spg/{A^{2 + sinh^o)), = u{t - to), (30) 
with asymptotic behavior at tq ^ oo 



lim|x(t)KA.exp(-|ro|), A, = 2{2y/'^su;/Ao. (31) 
On the basis of this solution one can prove that the quantity 

Xi = 1/(2^2) (32) 

and does not depends on / and t/. The final formula for det' {Lsad) entering in Eq. ( pOD for 



the partition function Z^pr can be easily obtained if we combine Eqs. 



det'iLsad) = det(Lo)/(2A2)". (33) 



(5). The quantity K from Eq. (|20D for the partition function Z^pj., which gives the 
measure of integration over time 0, can be calculated if we take into account the explicit 
form of the measure dGt in the functional integral (|2^ ) for the partition function for s ^ 1 



dGt = d'ilJt/V2^ ^ Kdt with K = sjAo/27r. (34) 
At present all quantities entering in Eq. (|20|) for Z^pr are defined and we have 

Z:^^ = (det(Lo))-V2(jT)"exp(-z$)/n!, (35) 



/ = Cexp(-Ao),C = V2AArK = 2 2^'^uoJs/ti 



where $ is the total Berry phase during all time of transitions. 
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For obtaning the final answer for Z it is necessary to sum all saddle point contributions 
to the partition function. These total saddle point contributions can be presented in the 
form of the following sum: 

Z = Zo[1 + (A + 5)(A + 5)TV2+ ... (36) 
[(A + fi)(A + 5)r2]7(2n!)+ ...]. 

where is the contribution of the trivial saddle point, A,B,A,B are instanton and anti- 
instanton without the determinant of the trivial saddle point. The dependence on the Berry 
phase can be easily determined: 

A, B = exp{iTTs)I , B, A = exp{—ins)I, (37) 
{A + B){A + B)= 4/2 cos^{tts). 

After these calculations we can sum all terms of the expression ( p^ ) in the form of 
cosh(— TA£'/2) and write the final answer for the energy of the ground state and its splitting 
by taking into account the explicit form of Zq (^S|) for the trivial saddle point: 

Eo ± 6E/2, Eo = Emin + (38) 
8 2^/^cjy^exp(-2sln(l + y2))cos(7rs)], 

--V2sg, g = f{2s-l)/s. 
We can see that this splitting is exponentially small for large s and equals zero in the 
case of half-integer spins in full agreement with Kramers theorem ||^J^ . The cancellation of 
the splitting for the case of the half-integer spin takes place due to the compensation of the 
contributions of the instantons of A and B type to the partition function Z. 

But, unfortunately, expressions ( p8D for the ground state energy Eq and the instanton 
splitting 6E are valid only qualitatively. The correct answer for Eq contains the correction 
—Sgs/A and for the instanton splitting 5E the additional factor 1/(1 + ^/2). The origin of 
these corrections is in the subtle nature of the functional integral which demands accurate 
treatments when we calculate functional determinants. 
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E± = 
5E = 

where Emin = u = 



IV. ENERGY OF THE GROUND STATE OR THE TRIVIAL SADDLE POINT 
CONTRIBUTION TO THE PARTITION FUNCTION 



The expression Eq for the ground state energy is wrong due to an incorrect definition of 
the product of operators at the same point in terms of the variables (p and p. The natural 
variables for the functional integral (§) are the z-variables. The trivial saddle points in terms 
of these variables are z = ±1. In the neighborhood of the saddle points we can represent 
the variables in the form 

z = ±l + ^2fs a, z* = ±l + ^J2fs a*. (39) 

The quadratic part of the action close to the saddle point is 

N-l 

^g(^) = [('^n+l'^n - Ki(^n) " " (40) 

n=0 

(/3/2)«2+, + a^)], 

where a = 3Tsg/N, (3 = Tsg/N. The quadratic form Ag{a) can be partly diagonalized if 
we pass to the cu-representation for the variables a„, a*: 

Af-l 

an = N~^^^ exp{iumn)am, (41) 

m=0 
N-l 



m=0 

where uJm = l-nimjN . The action in the tu-representation has the form: 

= XI [((1 ~ exp(-zt^„) - V)a*^ara - (42) 

m=0 



The Gaussian integrals with the quadratic form (|4l|) can be easily performed and the deter- 
minant can be obtained. The determinant D is equal to the product of the eigenvalues 



Am of the quadratic form ( ^T]) 



Am = + 2(1 -a)(l-coscu™), (43) 
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The partition function Z has the following expression in terms of Xm'- 



Z = exp(E„,„T)y, F = n A-^/'. (44) 

m=0 

For the calculation of the quantity Y it is convenient to represent the eigenvalues Am in the 
form 

Am = fimfJ'N-m, fim = b + Cexp^iu^rn), (45) 

(6, c) = (1/2) _ p2 ± ^(2 - ay - f3^^ 

- ±(1 - a/2 ± .\/a2 - /32/2). 

In terms of the quantities fj,m we have for F 

N-l N-1 

Y=l[fx-' = exp[- ln(/.m)] = (46) 

m=0 m=0 

exp[-iVln(6) - ^ ln(l + (c/6) exp(iia;m))]. 

The sum over m in Eq. (^) can be easily performed if we expand the logarithmic function 
in a series over x = (c/b) exp{iuJm.) and take into account only the terms with p = kN, 
where k is integer, which give nonzero contributions. The result for Y has the form 

Y = r^(l - {c/bf)-\ (47) 

Using the property a, 6, c ~ 1 + d/N we easily get for the partition function and the 
ground state energy 

Zo = 2exp[-EoT][l - exp(-u;r)]-\ (48) 

Eq = Emin + (l^ — to'o)/2, 

where uj = V2sg, cuq = 3sg/2. The factor 2 in this formula takes into account the presence 
of two minima in the action or the two types of the low-frequency excitations. 



We can see that the result (^) differs from the naive result Eq naive obtained previously. 



From the point of view of the functional integral this difference is due to the contribution 
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of the large eigenvalues of the quadratic form (1^) which correspond to the large uj in the 
partition function (P^. From the operator point of view the quadratic form (^) expressed 
in terms of operators (for example with the help of the Holstein - Primakoff representation 
with the quantization along x) has to be diagonalized with the help of the Bogoluibov's 
u — V transformation. As a result of such diagonalization the ground state energy is shifted 
by the amount {u — ujq)/2. 

V. CALCULATION OF THE DETERMINANT AND THE MEASURE OF 
INTEGRATION OVER TIME IN THE CONTINUUM APPROXIMATION 

A. Calculation of the determinant 



In this section we present the calculation of the ratio of the functional determinants 
(p!8|) Xi without the primary integration over the variable p which permits us to avoid the 
artificial limitation (p3D (^As^ ^ 1. Our consideration is based on the arguments of Coleman 
101 applied to the system of two differential equations of the first order in time. This system 
of equations can be easily obtained if we introduce the quantum corrections to the classical 

solution u° (P^) 

V't = v^^V^t, vt = y/s6ut. (49) 

The quantum correction to the action has the form 

rT/2 

A2 = -(1/2) / yuLl\,^at,dt, t,j = 1, 2, (50) 

J-T/2 

where ya = il^t, yt2 = Vt and the 2x2 matrix Lq can be presented in the form 

^o = -e|r + 5, e^ = -l, (51) 



e = 


(0 A 


, B = 


at 






^-1 0) 
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where the functions at, bt, Ct can be easily calculated if we take the variation of Eq. ([I5D 
with respect to ipt and uf 

at = gs{l + Ut) cos{2ipt), h = gsutsm{2ipt), (52) 
Ct = -gs{l + cos^{ipt)). 

For the trivial saddle point we have 

at = gs, bt = 0, Cj = -2gs. (53) 

One can see that the operator Lq is Hermitian. 



We can fix the class of eigenf unctions of the operator Lq (51) if we demand that 
?/'(— T/2) = iIj{T/2) = 0. For this class of eigenfunctions all eigenvalues of the operator 
Lq are nondegenerated and the ratio of the determinants Xi for the instanton saddle point 
and for the trivial saddle point is the ratio of the products of all eigenvalues. The same 
arguments as in leads us to the conclusion that 

Xi = C,*(2^/2)/C(2^/2)Ao. (54) 

Here ipinsti't) ^'trivi.'t) ^-^e eigenfunctions of the operator Lq (with zero eigenvalue) with 
the boundary conditions: 

V^°(-T/2) = and d^'^{-T/2)/dt = 1, (55) 

and Ao is the eigenvalue of the operator Lq which corresponds to zero mode. 

For the case of the trivial saddle point such solution can be constructed directly: 

y{t)^im,vit))= (56) 
u}'\smh{uj{t + T/2)), cosh(cj(t + T/2)/V2), 

where lu = \f2sg and 

V;L(2^/2) ^ exp(^T)/2^. (57) 
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For the case of the instanton saddle point the function 4'insti't) can be constructed on 
the basis of the zero mode eigenfuction of the operator Lq which is proportional to the time 
derivative of the instanton solution 



yt = V^sAi^l w?) = (58) 
2V2.A (- 2coshM V2sinh(ro) \ 

\ 1 + cosh (tq) cosh [To) J 

where A is the normalization constant due to the normalization condition 

/oo 
y^dt = 1. (59) 
-00 

This solution has the simple asymptotic form 

= 2ljA{-V2, ±1) exp{-\ro\) for |ro| ^ ±cx). (60) 

The solution (|58D does not satisfy the boundary conditions ( ^5|) and we have to construct 
the second linear independent eigenfunction of the operator Lq with the zero eigenvalue. 
The asymptotic form of this solution can be easily determined from the solution for the 
trivial saddle point 

y'^ = 2ujA{TV2, -l)exp{\To\) for |ro| ^ ±00. (61) 

It is important to note that these two solutions satisfy at any t the condition of conser- 
vation of the current or the Wronskian which can be easily deduced from the explicit form 
of the operator Lq 

I = yteyl = IGcuA^. (62) 
At this point we can construct the eigenfunction 4'insti't) 

ylAt) = -{AV2Aujr\e-^'^yt + e-^/^t). (63) 

One can easily check that ip^^^^iT /2) = 1. 

The last step which is necessary in order to determine the ratio Xi consists of finding the 
eigenvalue Aq which corresponds to the zero mode. For that it is convenient to convert the 
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equation on the eigenvalues A into the form of an integral equation. This integral equation 
has the form 

y\t) = yL.(t) + Ijy^iiy' ■ - (64) 

y1{yt'-y\t'))]dt'. 

Using conservation of the current / (|62D one can verify that this equation is equivalent to 
the original one. Moreover, the boundary condition %l)^{—T/2) = is satisfied. The second 
boundary condition in the framework of perturbation theory has the form 



IQu'^A^ J-T/2 

i^tiyt' ■ yLtimdt'. 



\ rT/2 

^ ^p\T/2) = 1 + Lm(y' ■ ylJt')) - (65) 



The integral in the right hand side of Eq.(^) can be calculated asymptotically 0] and we 
obtain the following relation for determination of Aq 

= 1- {\o/16V2uj^A'^uj) exp{ujT). (66) 



Combining Eqs. (|5^ , [57|j66|) we get for the ratio X; (|5 



Xi = ^-^ = (67) 

(exp(a;T)/2u;) ■ (A16y2cj3A2)exp(-cjT) ^ ^ 

1 



We take into account in Eq. (|67|) that the operator in the quadratic form ( pO|) is in fact 
proportional to A and all eigenvalues are also proportional to A. 

B. Calculation of the measure of integration over time 



The measure of integration over time can be determined if we follow basically [0]. The 
difference between our approach and ||^ is that we use the canonical formulation. The 
measure of integration over variables ip and v in the leading order (with respect to 1/s) is 
canonical 
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rf/i=n'^4^, AiV = r. (68) 



n=o v27rV27r' 

We stress in this expression for the measure that integration over each independent vari- 
able contains a factor 1/\/2ti. We want to pass to integration over the amphtude of the 
eigenfunctions of the operator Lq ([^). Then the measure takes the form 

2^-^ dc- 

d/i= n 4=^- (69) 



i=o v27r 



We want to find the connection between the amphtude of the zero mode cq and the time 
variable to which corresponds to the shift of position of the instanton. For that let us 
calculate the change of the original variables jt = (V^t? Vt) due to the change of Cq and to^ 



dyt = VAy^dco, (70) 
dyt = v^(v9°, u^t)dt = A-'y^dt. 

The additional factor -\/A preceding dco is due to the difference between continuos and 
discrete normalization. Comparing these two expressions we get 



dfio = dco/V2n = Kdt, K = A-^^s/2txA. (71) 

Therefore the constant C(^) which determines the instanton contribution to the partition 
function is 



C = X-'^^K = 22^l^^7fiiuj. (72) 

The normalization constant A completely disappears from the final expression for the con- 
stant C. This expression completely coincides with that one (|35|) which was obtained on the 
basis of paper for the action integrated out over the variable p. 

VI. CALCULATION OF CORRECTIONS TO THE DETERMINANT CAUSED 
BY THE DISCRETE NATURE OF THE FUNCTIONAL INTEGRAL OVER TIME 
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In our calculation we passed carelessly from the original variables z*, zt to the natural 
variables (pt, Pt- Small corrections of the order Au can lead to finite contributions to such 
quantities as the measure of integration and det(I/o). Practically they are the same quantities 
because by change of variables of integration they can be reduced to each other. At first we 
explain a mechanism of such corrections on a qualitative level and after that we apply it to 



our case. 



Let us suppose that each variable Cj in the measure (^) has a small multiplicative 
correction of the order Aa; 

^ dc- 

dfi^Y[^^{l + aiAu) =dfioR, (73) 



i=i V27r 

N 

R= Y[{l + aiAuj). 
1=1 

This factor R to dfio can be presented in an exponential form and it is different from unity 

R = exp{SAo), (74) 

N T 

deltaAo = y^\n{l + aiAuj) c::^ lu atdt, 
i=i 

where 5Aq is the correction to the classical action. If the function at —>■ const a.t t T then 
6Aq ~ const ujT . If the function a is localized over time within a region order uj~^ then 
5Aq ~ const . 

If we calculate the determinant the mechanism of influence of corrections is more subtle. 
Let us suppose that the matrix L can be presented in a form 

L = L^ + AujQ, (75) 

where the operator represents some discrete version of the operator Lq of (|5TD and the 
operator Q represents some discrete version of the second order differential operator 

Q = ao + aiuj + a2UJ ^ i — \ , (76) 

where functions oq, ai, 02 are some smooth functions of time t of the order of unity. Simple 
arguments shows that the contribution to the determinant due to the term is always 
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small contrary to the contributions of the ai and 02 terms which are of order of unity. To 
demonstrate this statement let us represent the operator L in the form 

L = [1 + AuQG^]L'', G° = (L°)-\ (77) 

and the determinant takes a form 

det[L] ~ det[L°] exp(-2Mo), (78) 
6A0 = -{1/2)AujTi[QG% 

One can check that in the discrete representation the matrix elements of the Green's 
function {G^)nn' are of the order A and the characteristic width over n — n' is of order 
(ujA)^^. This means that the contribution of the Oq terms of the operator Q into 6Aq is of 
order A. This can easily be checked if we convert the trace in Eq. (^) into an integral over 
t. 

The first impression is that the contributions of the ai and a2 terms to det[L] are also 
small because the Green's function {G^)nn' has a width over t — t' of order u;"^ and its 
derivatives have the order of magnitudes u. But these arguments are completely wrong 
because the Green's function is a singular function over t — t'. It has a jump at t = t' and 
the time derivative of this jump is of order A~^. If the functions ai and 02 —>■ const at t — > T 
then 6Aq ~ const ujT . If the functions ai and 02 are localized over time in a region of order 
oj'^ then 5Aq ~ const . 

At this point we can calculate corrections to the instanton contribution due to renormal- 
ization of the measure and the functional determinant in the leading order with respect to 
1/s. There are three sources of corrections of order A to the continuos approximation. The 
first origin of corrections is the Berry phase. The action with sufficient accuracy is of the 
form 



N-l 



1 + Z*^Zn (1 + Z*^ZnY 



(79) 



n=0 

The second origin of corrections is the Hamiltonian due to its dependence on 
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N—1 ( * *^/'Q _L * Q * 2 3\ 

Ahi = gs /X 2^ 7T^~m3 " ^^'^J 



The third origin of corrections is the difference between the average of the square of the 
Hamiltonian and the square of the of the average of the Hamiltonian: 



N-l 



= {i/2)A' y: 



n=0 







Zn) 


{Zn 


Zn) 



2- 



{Zn\Zn) 



(81) 



Because this correction depends (in the leading order) only on and z„ it can contribute 
only to the renormalization of the measure of integration. 

At this stage we can produce all calculations of the measure and the determinant for 
the instanton contribution with the necessary accuracy. For that we begin by interpreting 
the integration over the variables z'^ and z'^ {zn = z'^\ iz'^ and z'^ = z'^ — iz'^) in expression 
(^) for the partition function Z as an integral over the two-dimensional surface Im(z^) = 
and lm{z'^) = of the four dimensional space where variables and z'^ are considered as 
complex variables. Because the function (1 + z*z„)^^ exp(yl(z*, z„)) is an analytic function 
of the variables z'^ and z" it can be continued in the four-dimensional complex manifold and 
in this way we can arrive to the instanton saddle point. In the neighborhood of the saddle 
point we have to integrate over the two-dimensional manifold which realizes the directions 
of steep descent. The direction of steep descent is chosen correctly if all eigenvalues of the 
quadratic form in the exponent of the action are real. 

This program can be realized with the help of the following change of variables 2* and 
Zn in the neighborhood of the saddle point: 

Zn = Zn + {i/ y/s)Zn{'4'n " Vn/ (1 + uD), (82) 
Zn = Zn- ii/Vs)z*^{lpn + Vn/ (1 + M^))- 

Here we understand z„ and z* as classical (nonfluctuating) variables connected with the 
previously introduced variables (pn and by the relations: 



Zni^n, Un) = -y/(l - lUn) / {I + Wn) exp(iv5„), (83) 



zl{ipn,Un) = y^(l - iUn)/ (1 + m„) exp{-iipn), 
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thus the variables z„ and ^* are not in our case complex conjugated. We shall understand 
that the classical variables and or x„ = u„) satisfy the classical equations of 

motion which determine the saddle point with corrections of order A taken into account. 
This means that only in the main approximation the variables x are equal to x° = ((^9°, u^) 
determined by Eq. (IT3). The difference between variables x and variables x'' is of order A 



and is determined by small corrections of order of A contained in the action Ab (^) and 
corrections to the action Ahi ( ^Ol) and Ah2 (HD- It is fortunate that these corrections are 
nonessential to our problem. 

The reason for that is the canonical form of the measure (|68D in terms of the variables 
y„ = {ipni Vn)- Notice that the change of variables (0) was stimulated by formulas of 
differentiation over time of the quantities ^„((^„,m„) and ^*(v9„, 

At this stage we can calculate the renormalization of the functional determinant due to 
the A-corrections. This can be done on the basis of the following formula for the decompo- 
sition of the action in the neighborhood of the saddle point: 



A{z ,z)=Ao + - y ymVn'j + .. 



(84) 



ni,n'j ^3 

This formula strongly simplifies the calculations and can be proved if we consider the original 
variables Zn and z^ as function of variables yn = {'ipn, Vn) from one hand and the variables 
Xn = ifn, Un) from another. 

The matrix can be chosen on the basis in the form 

\ 



d+ + h„ 



-d- + h 

where the functions a„ , 6„ and c„ were determined in Eq. 
are determined by the following relations: 



(85) 

and the difference derivatives 



9_/„ = (/„-/„_i)/A, 9+/„ = (/„+i - /„)/A, 

d^fn = (/n+l + fn-l " 



(86) 



The Green's function satisfies the relation: 
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)n,naiG )na,n' — Sn,r, 



(87) 



and can not be found in general form. But this is unnecessary for our purposes. We are 
interested in the singular part of the Green's function (G°)n'n at n' ^ n. This singular part 
of the Green's function at n' ^ n can be found in a general form: 



/ 



Mn',n{n' - n)cn, On' ,n{^ - A(n' - 1 - n)h 

\^ -6'n'+i,n(l + A(ra' + 1 - Aen'+i,nin' - n)an J 

here 6n'n is the 6'-function defined in the following manner 



(88) 



o 

"n'n 



1, n' >n + l d+6n',n = A M„/_„ 



(89) 



0, n' <n d^On'+l^n = A ^Sn',n 

After some tedious calculations the singular (containing essential derivatives) part of the 
operator Q entering in Eq. ( [77| ) can be presented in a form: 

fn, (90) 



(Q/)„ = (1/2^) 



where is an arbitrary function. For the diagonal elements of the matrix Q it is sufficient 
to keep the second derivative only. For nondiagonal elements we have to keep the first and 
the second derivatives. Acting with the operator Q on the Green's function G'°)„',n and 
calculating the trace we get for the correction to the action 6Ao 

(91) 



5Ao = -(1/2) J {ft - Udt, 
ft = {l + u'iy^at - (1 + m2)q + A^tUt. 



Since we are interested in the ratio of determinants we subtract from the function ft its 
value at the trivial saddle-point /oo- Using Eq. (|l^) for (f we get the final expression for the 
correction to the instanton action 6Aq 



SAn 



-gs 



[1 - cos'^{ipt))dt = -ln(l + V2). 



(92) 



The obtained result is surprisingly simple. It can be found if we change in the expression 
for the energy splitting (^) s =^ s + 1/2. Such change has a quasiclassical meaning and can 
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be found by a simpler manner then in this section. Finaly,! result we have instead of the 
expression (pH]) for the ground state energy Eq and for the instanton splitting 

Eo = -s{s + 1 - V2)f, (93) 
6E = 16 23/^-1/^53/2/(1 + v^)-('^+i) cos(7rs). 

This result completely coinsides with the result of |^ for the case A = B = f . 

VII. COMPARISON WITH THE NUMERICAL RESULTS 

In this section we will compare the exact ground-state energy and splitting with the 
results obtained in the above discussion, namely eq(|38D. The Hamiltonian (|^) is easily 
diagonalized on the basis of the eigenfunctions of S"^ and 5*0. The ground-state belongs to 
the symmetric representation and, so, we will just consider this representation for different 
values of the total spin s. In table I the numerical results and the theoretical ones for the 
ground-state energy (^Sp and its splitting ( P^D are presented. As expected the relative error 
between the calculated and the exact ground-state energy decreases with s. 

VIII. DISCUSSION 

We want to discuss two points here: the interpretation of the effective change of s — »• 
s + 1/2 in the effective action and the procedure of calculation of corrections to the instanton 
approximation. 

A. Interpretation of renormalization of the instanton action 

Let consider the explicit form of the spin operators acting on the space of functions 
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-^mi^) = e^"''^/\/2^, < < 27r, m = -s, s. (94) 
The spin operators have the form [p|,p!0| 

S+ = ^{s + 1/2)2 _ _ 1/2)2 e^^^ 4=^^ (95) 
5_ = e-^^y^(s + 1/2)2 - - 1/2)2, p = -i9/9(/?. 

Substituting this representation for the spin operators (|95|) into the Hamihonian @ and de- 
composing it over p (this decomposition can be justified) we get in the leading approximation 
over 1/s: 

H{p, ^) = (//2)(p2(i + cos^ ^)-is + 1/2)2 ^^g2 (gg) 

We can calculate with this Hamiltonian the partition function applying the procedure of 
the p — ^ construction of the functional integral. In result we get the "functional integral" 
in which in each time section we have an integration over in the limits [0, 27r] and a 
summation over m in the limits [— s, s\. We can change summation over m on integration 
over p and prolong it to infinty. After integration over p we get the action which coincides 
with the action (|2^) with one essential difference that instead of in (^) we have (s + 1/2)2 
in (|96|). This means that we can get the correct answer for the tunneling splitting 5E ( |93D 
in the continuous representation 0]. Because we know at present that corrections are small 
(see discussion below) there is one question: why are corrections absent in the approach 
discussed in this section? The explanation lies in the difference between the coherent state 
construction of the functional integral and the p — ^p construction. One can check that at the 
calculation of the functional determinant in the p — ^ functional integral the corrections of 
the order A are absent. These corrections are also absent in the usual problem of tunneling 
in quantum mechnics which is confirmed by the coincidence of the result of the energy 
splitting with the quasiclassical one 0). 



B. Corrections to the leading approximation 
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We can construct perturbation expansion around of the trivial saddle point (^) as well as 
around of the instanton solution (0). Both these expansions are over very well defined small 
parameter 1/s. This parameter appears due to the presence of the chracteristic factor 
in the change of the variables (^) and ([8^). We want to stress that formulation in terms of 
the functional for the spin coherent states gives the regular procedure for the calculation of 
corrections. These corrections are of order const for the calculation the ground state enrgy 
and are of order 1/s for calculation of the factor before the exponent in the expression for 
the instanton splitting. 

In conclusion we want to stress one peculiarity of this perturbation theory: the presence 
in it of characteristic terms proportional of the number of the time sections A^. The presence 
of such terms is the characteristic feature of the perturbation theory when the measure of the 
integration is not trivial. In such theories the kinetic term is also nontrivial: the effective 
mass or the coefficient before a term with the time derivative is a function of the field 
variables. At the calculation in the framework of perturbation theory such terms lead to 
divergencies at large frequences uj. These divergencies have to compleatly canceal A^-terms 
which follow from the measure. 
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TABLES 

TABLE I. Results for the ground state energy and its splitting. In the first column is the 

magnitude of spin, the second and third cohimns are, respectively, the exact and the calculated 
ground-state energy, the fourth and the fifth columns give the exact and the calculated ground-state 
splitting. 
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